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Abstract 

Anomalous diffusion is discussed in the context of quantum Brownian motion with colored noise. 
It is shown that earlier results follow simply and directly from the fluctuation-dissipation theorem. 
The limits on the long-time dependence of anomalous diffusion are shown to be a consequence of 
the second law of thermodynamics. The special case of an electron interacting with the radiation 
field is discussed in detail. We apply our results to wave-packet spreading. 
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I. INTRODUCTION 



There has been a recurrence of interest in the phenomenon of anomalous diffusion. We 



refer to a samphng of the hterature 



erest m tnepnenomen 



There have also appeared in 



ihe literature a variety of explanations for such behavior. In particular, Tsallis and others 



a 



have proposed that conventional Gibbs-Boltzmann statistics be generalized by extending 
the definition of entropy. On the other hand, it has been shown that an explanation using 
conventional statistics is possible if one uses a Fokker-Planck equation with either fractional 
deHvative. Q or ti..depende„t diffusion coefficients 00. In addition, a stndy nsing path 
integral methods showed that anomalous diffusion appears in quantum Brownian motion 
with colored noise. 

Our purpose here is first of all to show that these same results follow directly from the well 
known exact formula for the correlation function that can be obtained using the fluctuation- 
dissipation theorem. In this section, we present general results which are applicable to an 
arbitrary heat bath. Then, in Sees. II to IV, we apply these results to the Ohmic heat 
bath, colored noise and a radiation heat bath (QED), respectively. In particularly, we note 
that the second law of thermodynamics places limits on the asymptotic behaviour. The 
physically realizable example of a charged particle coupled to the radiation field is given 
particular attention. In Sec. V, we apply our results to the spreading of a wave packet. Of 
special interest is the fact that the results in the QED case involve the bare mass. Finally, 
in Sec. VI, we present our conclusions. 

Diffusion of a Brownian particle is characterized by the long time behavior of the mean 
square displacement, which for motion during time t is defined to be 

s{t) = {[x{t) - x{o)r) . (1.1) 

Normal diffusion corresponds to a linear time-dependence for long times, with the diffusion 
constant given by = ^limt^oo s(t). Anomalous diffusion is characterized by a deviation 
from this linear time dependence of the form 

s(t)~t^+^, -1<7<1, (1.2) 



with 7 = being the case of normal diffusion. Our discussion will be based on the fami. 
formula for the correlation function, obtained using the fluctuation-dissipation theorem 



lar 
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1 h f°° fiuj 

- (a;(t)x(0) + x(0)a;(t)) = - / duj\m{a{uj + zO+)} coth — — cosujt. (1.3) 

2 TT JO 2kT 

Here a is the generalized susceptibility, which for a Brownian particle of mass m is of the 
form 

where fi{z) is the Fourier transform of the memory function. For our purposes all we need 
know about fl{z) is that, as a consequence of the second law of thermodynamics, it must be 
what is termed a positive real function That is, fl{z) must be analytic with positive 
real part everywhere in the upper half plane and, in addition, its boundary value on the real 
axis, which in general may be a distribution, must satisfy the reality condition, 

Re{fi{-uj + iO+)} = Re{fl{uj + zO+)} > 0. (1.5) 

Using the expression (jl.3|) for the correlation function, the mean square displacement can 
be written 

2h r°° Tioj 
s{t) = — dcjlm{a(c^ + i0+)}coth— — (1 -coswt). (1.6) 
TT Jo 2kT 

It is convenient to consider instead of s its time derivative 

sit) = — / duj\m{a{uj + iQ^)\uj coih. — — sincjt. (1.7) 
TT JO 2kT 

It should be clear from this formula that the long time behavior follows from the low fre- 
quency behavior of the integrand. We now show how this works for the most general case. 



and ()1.7|) becomes 



II. OHMIC CASE. 

Consider first the Ohmic case, where fi{z) = C,, a. constant. Then 

a(^ + zO+) = (2.1) 

2fi roo fujj 

s{t) = — rfw^-^— coth— — sinwt (2.2) 

The long time behavior is obtained by expanding the part of the integrand other than the 

sin cut function in powers of uo. That is, 

AkT , sin cut 2kT , , 

nC Jo u) C 
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Thus s is a constant and the mean square displacement is hnear in time. This is the case of 
so-called normal diffusion. 



2Ti r°° 2h 
s(t) ~ — : / dusmujt = -—. (2.4) 
vrC Jo nCt 



The case of zero temperature requires that in (|2.2|) we set coth 1 first and then 

expand in powers of u, to obtain 

2n ^ . 2h 
— duj sm ujt = —— . 
rC ^0 vrCt 

Note that this is a quantum result, vanishing in the limit h —>■ 0, while the finite temperature 
result ()2.3|1 is classical. This is a general feature: diffusion (normal or anomalous) is clas- 
sical at finite temperature, quantum at zero temperature. The fact that zero-temperature 
diffusion involves one power of t less than for finite temperature is also general. 

It will be of some interest to exhibit the exact zero-temperature result, which takes the 
form 



where 



POO ciri f 7/ ) x _ 

V(x) = / du \ ^ J = -[e-^'EUx) - e"Ei(-x)]. (2.6) 

JO + M"' 2 

From its form, it should be clear that the asymptotic form of the mean square displacement 
at zero temperature is given by 

.(t)~^(log^ + 7E), (2.7) 
where 7^; = 0.577- ■ ■ is Euler's constant. 

III. COLORED NOISE. 

Consider the case where in the neighborhood of the origin 

Il{z) = rnb^-^ {-izy , -1<7<1, (3.1) 

where 6 is a constant with the dimensions of frequency. It is easy to verify that this is a 
positive real function if and only if 7 is within the indicated range (we choose the branch 
with — vr < arg(2;) < vr ). With this, we see that (|1.7j) becomes 

.X 2h ^ cos(f7)coth|^ , , 

m = - duj . ,, sina;t. (3.2) 



TT 



m^[(f)i-7 + (A)i-7 + 2sin(f7)] 
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Expanding the integrand for small uo, we find for large t, 



^ ' Tcmb^-^ Jo cu'^+i 



Since 



we have the result 



1^ ~ 2r(l + 7)cos(f7)' ^ ^ ^ 

St — — -f<. 3.5 

m6i-^r(l + 7) ^ ^ 

Thus, we see that colored noise corresponds to anomalous diffusion. In addition we see that 
there is a physical reason for the restricted range of 7: it is a consequence of the second law 
of thermodynamics. 

The case of zero temperature follows in the same way as in the Ohmic case. Setting 
coth ^ — > 1 and expanding in powers of (j3.2|) becomes 

^ 2^a^ .00 ^^sin^ ^ -|^^^t-\ (3.6) 

Note that the case 7 = is not singular (lim^_^o cot(|7)/r(7) = 2/7r). Note also that this 
zero-temperature result involves one power of t less than the finite temperature result ()3.5j) . 

The extremes of the range of 7, 7 = ±1, require special attention. The case 7 = — 1 
corresponds to a generalized susceptibility of the form 

a{z) = ^ -. (3.7) 

This form corresponds to a free oscillator with force constant K = mP. That is, 7 = — 1 
corresponds to a harmonically bound particle. For this case one can obtain from ()1.7|1 the 
exact result 

, . h ^ hb ^ , , 

s(t) = — coth — — sm bt. (3.8) 
^ ^ m 2kT ^ ' 

Thus, with 7 = — 1 the particle is bound and there is no diffusion. This is consistent with 
(j3.5|) . for which the multiplier vanishes for 7 = — 1. 

The case 7 = +1 in p.ip clearly corresponds to a shift in the particle mass, with no 
dissipation. In the absence of dissipation, it should be obvious that at finite temperature 
the mean square displacement should grow with the square of the time, consistent with ()3.5|l 
for this case. There is, however, a non-trivial example for which dissipation is present and 
the mean square displacement at finite temperature nevertheless grows as the square of the 
time. This is the quantum electrodynamic (QED) case, which we discuss next. 



IV. QED CASE. 



In the case of a nonrelativistic electron coupled to the radiation field, the Fourier trans- 
form of the memory function in QED can be written in the form: IJ, |l5[ 

= 3^7T^' (^-^^ 
where f2 is a large cut-off frequency, related to the electron form-factor. With this, the form 
fll.4|) for the generalized susceptibility for a free particle becomes 

z + iQ 

where m is the bare mass and 

M = m + ^ (4.3) 

is the renormalized (physical) mass. We can use m as a parameter in place of fi, writing 

M — m 

Q = < m < M, (4.4) 

where 

Te = = 6.25 X 10-2^s. (4.5) 

With this, we can write the generalized susceptibility in the form 

M - m - iMzTe . , 

tt(^) = — , ^ o/ . ^ ■■ r- (4-6) 

Note that with this, m can be viewed as a measure of the strength of coupling, with < 
m < M. The limit m = M corresponds to no interaction (i.e., free particle) and the limit 
m = corresponds maximal coupling, where the cut-off has its largest value consistent with 
causality { Q = t^^) . 

We can rearrange the form ()4.6p to write 

1 (M -m)Te 

a(z) = 1 —- (A 7) 

^ ' Mz^ -iMz{M - m-imzTe)' ' 

Then we see that 

Im{a(^ + .0+)} = -^6'{u;) + ^, (4-8) 

M Muj[{M — my + [mujTe) \ 

where 5' is the derivative of the delta-function. With this, the expression fll.7|l for s becomes 

s{t) = 1 + duj- ^— ^ ^ -rcoth^smcjt. (4.9) 

^ ' M ttM Jo M - m 2 + mcuTe 2 2kT ^ ^ 
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A. Free particle. 



Consider first the case of a free particle, with no interaction with the radiation field. In 
this case only the first term in the expression (j4.9|) survives, and we have the result, 

kT 

s{t) = -jjf^, free particle, (4.10) 

valid for all times. This result follows from equipartition (<^|Mt>^^ = |A;T) which for a free 
particle is equally true in classical and quantum mechanics. Note that this is an exact result: 
for a free particle at zero temperature s{t) = for all time. 



B. Maximal coupling. 

Next consider the case of maximal coupling, where m = 0. In this case we find the exact 
closed form result [3], 

2kT 2hT, r°° , hcu 



sit) = 1 -\ / dti; coth — — sm a;t 

^ ^ M TiM Jo 2kT 



2kT , , vrm, 

Thus, for long time s grows linearly with time. Indeed the mean square displacement for 
long times is exactly of the form ()4.1()|1 for a free particle. This is what we should expect 
since in QED a particle moving with constant velocity feels no force. 
At zero temperature, the exact expression ()4.11|) becomes 

s{t) = — 4.12 

Note that for long times this result is an exception to the general rule that for colored noise 
the long time dependence at zero temperature involves one power of t less than at finite 
temperature. Here the difference is two powers of t. This phenomenon is already signalled 
in the general expression ()3.6|) for colored noise at zero temperature, where the multiplier 
vanishes for 7 = 1. 

We note that at short times the exact expression 1)4.111) diverges like t~^. This should 
be puzzling, since by it's very definition p.l|l s(0) = 0. This divergence is a quantum 
phenomenon, arising from the mass renormalization. In order to understand this, we consider 
next the general expression ()4.9p at zero temperature. 



C. General case at zero temperature. 



In the limit of zero temperature, the general expression ()4.9p becomes 

'^=^/ "^^ fM \2^f ^sm^t. 4.13 

ttM Jo [M — my + [muTe) 

For long times this vanishes like (at zero temperature the particle is at rest). For very 
short times, of order mrg/M, this increases rapidly from zero to a large value, after which it 
falls to zero again like t~^. In fact, this integral is exactly that appearing in the corresponding 
Ohmic case (|2.5p . We can therefore write 

. ^V((^l^). (4.14) 

where V{x) is given by (j2.6j) . This is an example in which the bare mass appears for 

extremely short times (i.e., extremely high frequencies) Another example is the canonical 

commutation rule Clearly, in the limit m — 0, ()4.14|) becomes exactly of the form 

fl4.12|) . but now with ^ interpreted as the principal value. As in the Ohmic case, in order 

to integrate to find the asymptotic form of the mean square displacement we need to know 

the behavior at small t. With the result ()4.14|1 we see that 

2hTe , (M -m)t , , 

sit) ~ — f log ^ ^ + 7i?. 4.15 

Thus, we see that the bare mass appears even in the long time behavior of s{t). We shall 
have more to say about this in our discussion of wave packet spreading. 



V. SPREADING OF A WAVE PACKET. 



As a somewhat different application of these results, we consider the spreading of a 
Gaussian wave packet. The situation is as follows. At time ti the Brownian particle, which 
before the measurement is in equilibrium at temperature T, is measured (i.e., detected) 
with a Gaussian instrument of width o"i and centered at Xi. Then at a later time ^2 a 
second measurement is made, again with a Gaussian instrument but now with width cr2 and 
centered at X2. Let W{1, 2)dxidx2 be the probability that the first measurement finds the 
particle in range dxi about xi and that the second measurement finds the particle in range 
dx2 about X2. The mean square width of the wave packet at time t2 is then given by 

/OO /'OO 
dxi dx2W{l,2){xi- X2y. (5.1) 
-OO J — OO 
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Here we must be sure to distinguish this expectation from that in here xi and X2 are 

instrumental parameters (i.e., c-numbers) not quantum mechanical operators. One could 
well argue that for a quantum particle this wave packet width is a better measure of 
anomalous diffusion than the mean square displacement s. The picture, to repeat, is the 
following. An initial wave packet of width ai is formed at time ti. In the course of time this 
wave packet will spread. It's center will not move since the mean velocity is zero. At a later 
time t2 the mean square width, including the width a2 of the second instrument, is given 
by (EH). 

Now, for quantum Brownian motion an exact expression for W{1, 2) can be obtained jl7l |. 

^(1'2) = ^_exp{-— ^(4 - 2^ + 4)}, (5.2) 

where (note a misprint in Eq. 7.18 of [l"^) 

= af + (x{tif 



^2 ^2,/^., x2\ Wi),x{t2)f 
T = (72 + \X[ti 



2 



(TTp = ^ {x{ti)x{t2) + X{t2)x{ti)) . (5.3) 



With this, we find that ()5.H) becomes 

w\t2 - h) = s{t2 - ti) + al - "y^^^^f^^"^' + ^2- (5.4) 

Thus, the mean square width of the wave packet is a sum of four terms. The first is the 
mean square displacement p.l|) . The second and third terms involve the finite width of 
the initial wave packet, with the third term clearly arising from the uncertainty principle. 
Finally, the last term is the mean square width of the second instrument. Clearly, this final 
width can be taken to be zero, but not the initial width. 

For quantum Brownian motion, the commutator is a c-number given by [l^ 

[x{ti),x{t2)] = / dujlm{a{uj + i0'^)}smuj{t2-ti). (5.5) 

TT JO 

Note that the commutator is temperature independent. We now consider a number of cases. 
A. Colored noise 

For long times, using the method described above, we find for colored noise, 



where to save writing we have written t = t2 — ti. At finite temperature, using this result 
for the commutator and the corresponding result ()3.5|) for the mean square displacement, 
we see that for — 1 < 7 < 1 and for sufficiently long times the term s{t2 — ti) dominates 
in (j5.4j) . Thus at finite temperature, the long time behavior for colored noise is the same, 
whether measured by the mean square displacement (ll.lj) or by the wave packet width (j5.1|) , 
and always classical. 

On the other hand, at zero temperature, using the result ()3.6|) we see that the general 
expression ()5.4|) becomes 

r(l + 7)m6i-T 

Here the first term corresponds to s{t) and we see that for — 1 < 7 < this term still 
dominates at zero temperature. The Ohmic case 7 = appears to be singular, but as we 
noted above this is not a problem. In fact, for the Ohmic case at zero temperature we have 



— log — + (Ti + 



w^OhmicW ~ — :log— + ^1 + 7T27^ + ^2- (5-8) 



Thus, in the Ohmic case, we see that the first term still dominates. On the other hand, for 
< 7 < 1, for sufficiently long times the third term dominates, corresponding to what we 
might call uncertainty principle spreading. 



B. Free particle 

For the free particle, using the expression (j4.8p with m = M, the expression (|5.5|) for the 
commutator becomes 

[xih),xih)]='^ih-h). (5.9) 

With this and the expression ()4.10|) for the mean square displacement, we see that the 
general expression ()5.4|1 becomes 

= ^«= + <^? + j|^ + -i. (5.10) 

First of all, we note that at zero temperature and with cr| = this becomes exactly the 
formula for wave packet spreading found in textbooks For finite temperature the 
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first term, the mean square displacement, has the same dependence as the third term, 
the uncertainty principle spreading. The ratio of these two terms is A^/Svrcrf, where A = 
J2TTh^/MkT is the mean thermal de Broglie wavelength. Thus, we see that the spreading 
will be classical provided that A <^ ai. 



C. QED 

As in our discussion of colored noise, the case 7 = 1 requires special consideration and 
we discuss it in terms of the QED case. For the general form ()4.8|) of Im{a}, we see that 
()5.5|) becomes 

[x{h),x{h)]QED = jT{t+- duj^^ ^ — }. 5.11 

M TT Jo u![{M — my + {7nujTf,y\ 



The integral appearing here is standard and we find 

[x(ti),x(t2)]QED = ^{t + re(l-e-(^-'")*/'"-)}. (5.12) 

Here the exponential is exceedingly small for any time of the order of or longer. It does, 
however, guarantee that the canonical equal-time commutator, [a;(t), x(t)]QED = ih./m, does 
involve the bare mass, as it must [lfi| . 

With this result, it should be clear that for long times the leading behavior of the wave 
packet width in QED is the same as that for a free particle. There will be lower order 
corrections, of course, but these are perhaps best discussed in terms of the zero temperature 
result. At zero temperature, using the results ()4.15|) and ()5.12|) we find 

VI. CONCLUSIONS 

The first calculation of normal diffusion [equivalent to our equation (2.3}1 was carried out 
by Einstein in his famous explanation [2^ for the observations of Brown j2l| on the random 
motion of pollen grains immersed in a fluid. The term "Brownian motion" is now used 
in a generic sense to denote random motion and it covers a wide spectrum of phenomena. 
Einstein's explanation of the original Brownian motion used a diffusion equation but, shortly 



afterwards, Langevin 



22| | presented an entirely new approach to the subject, based on the 
11 



introduction of a stochastic differential equation, wliich in the words of Chandrasekhar , 
constitutes the "modern" approach to this and other such problems. Langevin's analysis 
was confined to the classical long-time motion of a free particle at high temperature in an 
Ohmic heat bath (based on Stokes analysis of viscous drag). 

Here, we use a generalized Langevin equation which incorporates quantum and non- 
Markovian (memory) effects, arbitrary times, arbitrary temperature and a very general heat 
bath. This leads to a fluctuation-dissipation theorem which enables us to calculate anoma- 
lous (in a sense that it goes beyond the Einstein-Langevin case) diffusion for a wide variety 
of parameters which, in turn, will facilitate analysis of a variety of physical phenomena. 
In particular, we applied our results to wave-packet spreading. We consider the use of the 

simpler and more 



quantum Langevin equation and the fluctuation-dissipation theorem to be 
physically appealing than approaches based on path-integral techniques j9|,|2J]. In addition, 
our description of the motion of a quantum particle in a linear passive heat bath is char- 
acterized by Re[p,{z)], the Fourier transform (spectral distribution) of the memory function 
appearing in the quantum Langevin equation. As a consequence of the second law of ther- 
modynamics, fL{z) must be a positive real function Q], that is fi{z) must be analytic with 
positive real part everywhere in the upper half plane and, in addition, its boundary value 
on the real axis, which in general may be a distribution, must satisfy the reality condition, 
given by (1.5). As a consequence, this limits the forms of fi{z), as exemplified in (3.1), which 
in turn limits the forms of anomalous diffusion. Finally, we discussed in detail the special 
case of an electron interacting with the radiation field. 

Acknowledgments 

We wish to thank the School of Theoretical Physics, Dublin Institute for Advanced Stud- 
ies, for their hospitality. 



[1] J. Klafter and I.M. Sokolov, Phys. World Aug. 2005, p. 29. 

[2] M.F. de Andrade et al., Phys. Lett. A 347, 160 (2005). 

[3] E. Lutz, Phys. Rev. E 64, 051106 (2001). 

[4] C. Tsallis and D.J. Bukman, Phys. Rev. E 54, R2197 (1996). 

12 



[5] M. O. Scully, G. Sussmann and C. Benkert, Phys. Rev. Lett. 60, 1014 (1988) 
[6] M. O. Scully, M. S. Zubairy and K. Wodkiewicz, Optics Comm. 65, 440 (1988) 
[7] C. Benkert, M.O. Scully, A.A. Rangwala and W. Schleich, Phys. Rev. A 43, 1503 (1990). 
[8] S. Schaufler, W.P. Schleich and V.P. Yakovlev, Phys. Rev. Lett. 83, 3162 (1999). 
[9] H. Grabert, P. Schramm and G.-L. Ingold, Phys. Rev. Lett. 58, 1285 (1987). 
[10] H.B. Callen and T.A. Welton, Phys. Rev. 83, 34 (1951). 

[11] G.W. Ford, J.T. Lewis and R.F. O'Connell, Ann. Phys. (NY) 185, 270 (1988). 

[12] G.W. Ford, J.T. Lewis and R.F. O'Connell, Phys. Rev. A 37, 4419 (1988). 

[13] Tables of integral transforms, edited by Staff of the Bateman Manuscript Project, Bateman 

manuscript project. Vol. 1 (McGraw-Hill, New York, 1954). 
[14] G.W. Ford, J.T. Lewis and R.F. O'Connell, Phys. Rev. Lett. 55, 2273 (1985). 
[15] G.W. Ford and R.F. O'Conneh, Phys. Rev. A 57, 3112 (1998). 
[16] G.W. Ford and R.F. O'Conneh, Journ. Stat. Phys. 57 (1989) 803. 

[17] G. W. Ford and J. T. Lewis, Probability, Statistical Mechanics, and Number Theory, Advances 

in Mathematics Supplemental Studies, Volume 9 (Academic Press, Orlando, Florida 1986) 
[18] G.W. Ford, J.T. Lewis and R.F. O'Connell, Phys. Rev. A 64, 032101 (2001). 
[19] L. L SchifF, Quantum Mechanics (McGraw-Hill, New York 1949) p. 58. 

[20] A. Einstein, Eine neue Bestimmung der Molekilldimensionen, (Wyss, Bern, 1905) [Einstein's 
thesis]; ibid. Ann. d. Physik 17, 549 (1905); this and the further articles of Einstein appear 
in Investigations on the Theory of the Brownian Movement (Dover, New York, 1956). 

[21] R. Brown, Philos. Mag. 4, 161 (1828); Ibid. 6, 161 (1829). 

[22] M. P. Langevin, C. R. Acad, Sci., Paris 146, 530 (1908) (a translation of this article appears 

in D.S. Lemons, A. Gythiel, Am. J. Phys. 65 1079 (1997)). 
[23] S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943). 

[24] V. Hakim and V. Ambegaokar, Phys. Rev. A 32, 423 (1985). So far as we know, this paper 
contains the first derivation of the width of a spreading wave packet, equivalent with our Eq. 
(5.4), but only for the ohmic case. Also, we note a misprint in their expression (38) for the 
width. 



13 



